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ABSTRACT 

Galaxy clusters will distort the pattern of temperature anisotropics in the microwave background 
via gravitational lensing. We create lensed microwave background maps using clusters drawn from 
numerical cosmological simulations. A distinctive dipole-like temperature fluctuation pattern is formed 
aligned with the underlying microwave temperature gradient. For a massive cluster, the characteristic 
angular size of the temperature distortion is a few arcminutes and the characteristic amplitude a few 
micro-Kelvin. We demonstrate a simple technique for estimating the lensing deflection induced by 
the cluster; microwave background lensing measurements have the potential to determine the mass 
distribution for some clusters with good accuracy on angular scales up to a few arcminutes. Future 
high-resolution and high-sensitivity microwave background maps will have the capability to detect 
lensing by clusters; we discuss various systematic limitations on probing cluster masses using this 
technique. 

Subject headings: cosmology: theory - galaxies: clusters: general - cosmic microwave background 
gravitational lensing 



1. INTRODUCTION 

In the post-WMAP era of microwave background mea- 
surements, attention is quickly shifting to smaller an- 
gular scales. At scales below 4 arcminutes, the tem- 
perature fluctuations are dominated not by primordial 
fluctuations associated with the last scattering surface 
but rather by secondary fluctuations induced by inter- 
actions with the matter distribution at lower redshifts. 
Most prominent among these effects are th e Sunyaev- 
Zeldovich effect ijSunvaev &; Zerdovichlll972]) and gravi- 
tational lensing. 

This paper studies the gravitational lensing of the 
microwave background by galaxy clusters. Much re- 
cent effort has been devoted to lensing of the mi- 
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radio sources by galaxies, galaxy clusters, and large-scale 
structure is a major theoretical and experimental indus- 
try. Curiously little effort so far has been put into 
lensing of the microwave background by clusters; the 
onl y major work on the subje ct is the seminal paper 
by Efiliak & Zaldarriagal §000). The likely reason for 
this benign neglect is that only recently have microwave 
maps of sufficient resolution and sensitivity to detect 
cluster lensing becom e a realistic expectation l|Kosowskvl 
120031 iWoottenl 12002ft. The subjec t has recently been 
taken up b v iDodelson fc Starkmanl l)2003j) as it relates 
to galaxies; iCooravl l)2003|) explicitly showed the depen- 
den ce of the len s ing si gnal on cosmological parameters, 
and lBartelmannl (j2Q03) showed a map of a cluster lensing 
the CMB as an application of numerical techniques for 
gravitational lensing. 

Observing cluster lensing of the cosmic microwave 
background requires temperature variations on scales 



large compared to the cluster (so that we have some- 
thing to lens), but minimal variations on scales compa- 
rable to the size of the cluster (so that the lensing signal 
can be cleanly separated from the intrinsic temperature 
fluctuations). Notably this is exactly what the primary 
microwave background fluctuations offer. A microwave 
background temperature map is dominated by fluctua- 
tions on scales of a degree or larger, arising from density 
and temperature perturbations at the surface of last scat- 
tering at a redshift z ~ 1100. At scales smaller than a 
degree, the power spectrum of temperature fluctuations 
begins to decline due to diffusion damping: perturbations 
on scales smaller than the thickness of the last scattering 
surface are exponentially suppressed. On typical cluster 
angular scales of a few arcminutes, the primordial tem- 
perature perturbations are negligible. Secondary, non- 
linear temperature fluctuations will arise, but these are 
generally either tiny (e.g., gravitational lensing by large- 
scale structure) or have a frequency dependence differ- 
ent from the blackbody distortions due to lensing (e.g. 
the thermal Sunyaev-Zeldovich effect). The only source 
of significant nearly-blackbody fluctuations on scales of 
galaxy clusters besides lensing is the kinematic Sunyaev- 
Zeldovich distortion, mainly due to the peculiar motion 
of the cluster itself as well as a sm all component from 
large scale structure l|Vishniacll987|) . In general, the kSZ 
signal has a different morphology than the lensing signal, 
and in most cases the difference between the two will be 
clear. Furthermore, most of the kSZ signal should be 
strongly correlated with the thermal SZ effect. However, 
the kSZ signal will be a large source of contamination in 
the central regions, with a significant source of noise com- 
ing from bulk motions within the cluster due to objects 
that have recently been accreted. 

Seljak and Zaldarriaga (2000) considered idealized 
spherical clusters lensing a pure temperature gradient. 
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They also provided an extensive list of potential system- 
atic effects which must be overcome to observe the signal. 
In the space of three years, the sensitivity and resolution 
of envisioned microwave background measurements have 
increased dramatically, prompting a more detailed as- 
sessment of the cluster lensing signal. In this paper, we 
use model mass distributions from clusters in a cosmo- 
logical N-body simulation to lens a background Gaus- 
sian temperature field constructed from the temperature 
power spectrum in a realistic cosmological model. The 
cluster lensing signal is clear in the resulting map. Sec- 
tion II displays how the signal varies with cluster mass 
and cluster location on the sky, and Section III considers 
how well the cluster lensing signal can be inferred from 
an ideal map using a correlation function technique for 
estimating the unlensed temperature distribution. More 
elaborate model-fitting techniques are likely straightfor- 
ward but beyond the scope of the paper, requiring ex- 
plicit and careful considerations of specific instruments 
and observing strategies. Finally, we discuss potential 
advantages and systematic limitations of this method for 
cluster mass determination in the context of realistic ex- 
periments. 

2. ORDER- OF-MAGNITUDE ESTIMATES 

The root-mean-square temperature gradient in the pri- 
mary microwave background fluctuations is around 12 
/xK per arcminute for standard cosmological models con- 
sistent with the measured microwave power spectrum; 
the magnitude of the lensing signal is on the order of 
the local temperature gradient times the characteristic 
angular deflection. Typical deflection angles (roughly 
the cluster gravitational potential in units of c 2 , $ /c 2 ~ 
10~ 4 ) are on the order of tens of arcseconds. The frac- 
tional perturbation due to lensing is therefore on the or- 
der of a few fj,K. 

This is at least an order of magnitude larger than 
pert urbations due to the tran sverse velocity of the clus- 
ter ifBirkinsha w fc Gulilll983|) . which is on the order of 
the cluster potential <i>/c 2 times the transverse velocity 
vtrans/c. Typical peculiar velocities are a few hundred 
km/s, so this term should be at the level of roughly 
0.1 /xK. Similarly, anisotropy imprinted from the clus- 
ter potential changing over the photon travel time (the 
Integrated Sachs- Wolfe effect in the non-linear regime) 
should be roughly given by the potential times the pho- 
ton travel time relative to the dynamical time. With 
travel times on the order of 1 Myr and dynamical times 
on the order of 1 Gyr this term should be compara- 
ble to the transverse velocity effects but well below the 
expected level of lensing effects. Therefore, the signa- 
tures of moving or dynamically evolving clusters will be 
strongly suppressed compared to lensing of the CMB, 
even if the thermal and kinematic SZ effects can be effi- 
ciently removed. 

As discussed below, the kinematic SZ effect is also ex- 
pected to be on the order of a few /iK, with a similar 
spectral dependence to gravitational lensing. These two 
signals must be separated via their spatial morphology 
or other techniques. 

3. LENSED MAPS 

Figure ^ displays the lensing of a pure temperature 
gradient by a numerically simulated cluster of galaxies of 
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Fig. 1. — Gravitational lensing of a pure background gradient of 
15 fj,K per arcminute by a galaxy cluster of mass 7 X 10 14 h~ 1 MQ 
extracted from a numerical simulation. Background contours are 
spaced by 5 fiK and show the unfiltered lensed signal, while solid 
contours, spaced by 0.5 fiK, show maps that have been high-pass 
filtered, with all k < 1150 (i < 7200) removed. Each image is 5' x5' 
and does not include kinematic SZ effects. Blue (thick) contours 
indicate positive temperature differences, while red (thin) contours 
indicate negative values. 



mass M = 6/i _1 x 10 14 M Q at a redshift of z = 0.5, while 
Fig. shows the same cluster lensing two realizations 
of a Gaussian random temperature field with the power 
spectrum of the microwave background for a flat ACDM 
cosmological model with n = 1, h = 0.7, fl^h 2 = 0.024, 
and f^A = 0.7. Equal-temperature contours are plotted 
at 5 /iK separations as light dotted lines. To display 
the lensing signal more clearly, we then high-pass filter 
the map with a 3-arcminute filter scale. The tempera- 
ture contours of the filtered map are plotted at 0.5 [iK 
separations as heavy lines. 

The cluster i s dra wn from the VIRGO 1 simulations 
l)da Silva et alJ 1200(1) . Outputs (including gas) were 
available at z = 0, and for our lensing calculation we 
artificially placed the simulation volume at z — 0.5. To 
be conservative, we did not scale the box size to account 
for the expansion of the universe between z = 0.5 and 
z = 0. This scaling would lead to an overestimate of 
the central concentration of virialized objects at z = 0.5 
in the simulation, and therefore an overestimate of the 
lensing signature. The box was translated such that a 
massive cluster was at the center of the projected mass 
distribution, and the dark matter distribution was pro- 
jected onto a 2048 x 2048 grid using a simple nearest grid 
point method. The gas was smoothed (in projection) 
over the mean distance to the 24th nearest neighbor as- 
suming a uniform disk smoothing kernel. The resulting 
surface mass density map was used to generate deflec- 
tion maps with resolution of roughly 10" , using an FFT 

1 http://virgo.sussex.ac.uk 
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Fig. 2. — Same as Fig. 1, except that the cluster is lcnsing two 
different realizations of a CMB sky. 



method. 

We compute the reduced angular deflection vector due 
to the cluster mass distribution in the thin-lens approx- 
imation, 

a(x) = i / K (x')^^dx' (1) 

7T J |x — x'| 

where x is a 2-dimensional vector describing the angular 
position on the sky, k(x) is the projected convergence 
related to the surface mass density £ by 



k(x) 



AttG D d D d , 
c 2 £>„ 



E(x) 



(2) 



with Z?d, D s , and D ds the angular diameter distances 
from the observer to the lens, the observer to the source, 
and the lens to the source, respectively. The actual 
deflection vector a is related to the reduced deflection 
by the distance scali ng a — (Dd s /Ds)ct. (See, e.g., 
iSchneider et alJ l|1992f) for a detailed exposition of stan- 
dard gravitational lensing theory.) 

A lensed map can be constructed from an unlensed 
background map and a surface mass distribution E(x) 
by solving the lens equation 



/3(x) 



a(x) 



(3) 



where (3 is the source position in the image plane. For a 
source covering the entire sky like the microwave back- 
ground, the lens equation is straightforward to solve by 
pixellating the image plane, then using the displacement 
vector to map each point in the image plane to a point of 
the source plane, then assign the source temperature at 
that point to the original point in the image plane. The 
resulting lensed maps in Fig. ^ represent the signal an 
ideal experiment with very high angular resolution and 
sensitivity might measure. 

The CMB maps were generated using a CMB power 
spect rum generated by CMBFast 2 ( Seliak & Zal darria"gal 
f996). To be conservative, we use the lensed power spec- 
trum, which has slightly more power on small scales, but 
assume the underlying CMB anisotropies are a Gaussian 
random held. We construct unlensed maps over small 
regions of sky using the flat sky approximation. 

The cluster lensing signals are obvious in the filtered 
maps. It is clear that lensing induces structure in the 
maps on significantly smaller scales than those on which 
the microwave background has significant temperature 
fluctuations. For a cluster situated in front of a region 

2 www.cmbfast.org 



of the microwave background which is approximately a 
temperature gradient, lensing produces a characteristic 
dipole-like pattern, with a cool and a hot peak. The 
peak-to-peak amplitude is proportional to the gradient 
magnitude, but is generally on the order of 1 to 10 fiK; 
the angular separation of the peaks is on the order of 
an arcminute, with noticeable lensing effects out to radii 
of several arcminutes. Also note that a vector from the 
hot to the cool lensing peak must be in the direction of 
the local CMB gradient: this distinctive signature can be 
used to discriminate between a lensing signal and other 
effects local to the galaxy cluster, which will be physically 
uncorrelated with the microwave background. 

The gross features of the cluster lensing signal depend 
on cluster mass, redshift, and sky location. The lensing 
effect depends on the angular diameter distance to the 
cluster; for massive clusters which are distant enough 
that their angular size is on the order of an arcminute, 
their angular diameter distance is only a weak function 
of redshift. The lensing deflection depends primarily on 
the total mass of the cluster. The other crucial factor 
in the cluster lensing pattern is sky position: the am- 
plitude of the lensing signal is proportional to the local 
background temperature gradient. To extract mass in- 
formation about the cluster, we must be able to infer the 
unlensed background temperature pattern. This is eas- 
iest when the cluster is situated in a sky region which 
has a fairly uniform temperature gradient; in this case 
the characteristic dipole lensing pattern is produced. If 
the cluster is located in a region of the sky where the 
temperature field is near an extremum, or where the 
isotemperature contours are strongly curved, extraction 
of the lensing mass is more uncertain, because the un- 
lensed temperature field is not as well constrained by the 
lensed map. 

Figure|3shows clusters with a range of different masses 
lensing the microwave background. Here the simulated 
lensing mass distribution is a cubic volume 100/i -1 Mpc 
on a side drawn from the VIRGO simulations, placed 
at a redshift of z — 0.5. Both the amplitude of the 
lensing peaks and their angular separation increase with 
cluster mass. To gauge the effects of sky position, we 
have outlined with a square the location of all clusters 
in the simulated mass distribution with M > 10 14 M Q . 
Some marked clusters do not obviously correspond to 
clean features in the filtered map, indicating that the 
lensing signature is very difficult to observe at these po- 
sitions. Roughly half of the I0 14 M clusters in Fig. [31 
are in sky positions leading to appreciable lensing signals. 
As a rule of thumb, the lensing signal will be easiest to 
estimate when |VT| is large and V 2 T is small. At the 
same time, occasional less massive objects that are par- 
ticularly well-placed at the positions of strong gradients 
display clear lensing signatures, as with the cluster near 
(—30', 5'). Note that the positions of all clusters will be 
known from their thermal Sunyaev-Zeldovich distortions, 
which are at least an order of magnitude larger than the 
lensing signature. 

In Figure 0] we show the problems that can be intro- 
duced by significant substructure in the kSZ signal. This 
case is chosen as an unfavorable kSZ signal, although it 
could be even less favorable if aligned exactly with the 
gradient. Some of the kSZ is removed by selecting only 
that component that is anti-symmetric around the clus- 
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Fig. 3. — 50' X 50' field showing several clusters lensing the CMB. 
All positions of clusters with masses above 10 14 /i — 1 Mq are shown 
as open squares in the map, and background dotted contours are 
spaced every 15 fiK. Solid contours show a high-pass filtered map 
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Fig. 4. — Effects of kinetic SZ signal: left two panels show 
case of no kSZ, right two panels include a large contaminating kSZ 
signal. Top panels show a simple high-pass filter, as in figure 1, 
while bottom two panels show the effect of a high-pass filter that 
has also selected the component that is anti-symmetric around the 
cluster center (i.e., the imaginary part of the Fourier transform). 



ter center (as the lensing signal should be). Along the 
direction of the gradient, it is still possible to see the 
characteristic lensing distortion, but it is clear that the 
central region will be useless for the purposes of lens- 
ing reconstructions. We emphasize that this was not 



selected as a typical contamination, and the simulation 
physics (no feedback or preheating) led to gas cores that 
are much more concentrated than in observed clusters. 

4. DEFLECTION ESTIMATION 

The maps in the previous section make clear that in 
many cases (but not all), some estimate of cluster mass 
within a certain radius can be made from the pattern of 
lensing in the microwave background. It is likely that 
some kind of model fitting for both the cluster mass pro- 
file and the unlensed temperature map will give the most 
accurate cluster mass recovery, but a complete paramet- 
ric study is beyond the scope of this paper. Here we 
present a crude technique for estimating the lensing de- 
flection in the inner region of the cluster which works 
well in many cases. This deflection angle can then be 
used as a constraint on the mass distribution. 

4.1. Wiener Estimation of Unlensed Background 

For the case of a source intensity distribution which 
is a pure gradient, it is straightforward to estimate the 
lensing deflection in the direction along the temperature 
gradient. At sufficiently large distances from the clus- 
ter center, the background gradient will be clear and can 
be robustly estimated. For the case of the microwave 
background, a pure gradient is generally not a good de- 
scription of the background temperature distribution on 
scales of a few arcminutes over which lensing deflections 
are significant. However, the background temperature 
distribution is close to a gaussian random field, and the 
statistical nature of this distribution can be used to esti- 
mate the unlensed background temperature distribution 
in the region of a cluster: we mask out the central re- 
gions of the image, where lensing is known to be impor- 
tant, and interpolate an approximate unlensed temper- 
ature distribution using the region that is more than a 
few arcminutes away from the center. 

The most reliable method we found is Wiener interpo- 
lation. We assume that the correlation function of the 
microwave background temperature is known from power 
spectrum measurements. Then we have a good estimate 
of (TiTj) — C(\9i — 9j\) = Cij, the mean product of tem- 
peratures at two points 9i and 6j in the map; here C{9) is 
the temperature correlation function for two points with 
angular separation 9. The matrix C is just the theo- 
retical gaussian random correlation matrix constructed 
from the power spectrum. Given accurate estimates of 
the temperature map outside the cluster region, we can 
use these correlation functions to interpolate across the 
cluster. 

The pixels of some observed map can be viewed as a 
single data vector d. We then form the noise-weighted 
covariance matrix D: this covariance matrix includes the 
theoretical covariance C plus the covariance arising from 
pixel noise. We mask out given pixels in the region of 
a cluster by setting the diagonal elements of D corre- 
sponding to those pixels to a value which is large com- 
pared to the other elements of D. We then estimate the 
temperature map t es t defined by the theoretical correla- 
tion function which matches the measured pixels in the 
region ouside the cluster as t es t = CD M. This con- 
struction requires inversion of a large matrix. For a map 
with N x N pixels, the computation of the inverse co- 
variance matrix scales as A 6 , but pixels well beyond the 
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Fig. 5. — Reconstruction of unlensed CMB background field 
using Wiener interpolation. The axes are in arcseconds, and the 
diagonal line is marked in units of 5 arcseconds; temperature con- 
tours are spaced at 5 [xK intervals. Reconstruction of the central 
5' X 5' section of the map is done using the displayed map region 
with an effective pixel scale of 10" . Solid lines show the input map, 
dashed lines show the reconstructed map. 



cluster region do not contribute much to the cluster re- 
gion. Also, note that this inversion must be done only 
once for a given microwave background power spectrum, 
pixel scale, and mask geometry. 

As a verification of the method, we start with a CMB 
realization (with no cluster lensing), cut out the central 
8' and use the Wiener method to reconstruct the missing 
pixels. Results are shown in Figure[S] The reconstruction 
is not perfect, but the reconstruction errors are on the 
order of 5" , small compared to typical lensing deflection 
angles of 30-40" . For Fig. [5] we used a map that was 
64 x 64 pixels with a 10" pixel scale. This reconstruction 
took a few minutes on a desktop computer without any 
particular optimization effort. 

We also experimented with fitting a simple gradient or 
two-dimensional splines of varying stiffnesses, with vary- 
ing degrees of success on a case-by-case basis. For the 
Wiener interpolation technique, assuming the outer re- 
gions are unaffected by lensing introduces an unnatural 
boundary condition, so an iterative approach would im- 
prove this estimate. Ultimately, a simultaneous fit for 
the unlensed CMB distribution, kinematic SZ, and de- 
flection angle will be required for the most reliable con- 
straints, and will use the corresponding thermal SZ map 
to provide a rough template for the kinematic SZ signal. 

We note that the Wiener estimation technique pro- 
vides a general method for testing Gaussianity of the 
microwave background temperature distribution. 

4.2. Deflection Angle Reconstruction 

Applying the Wiener estimation method, we show the 
reconstructed deflection angle along a slice through the 
cluster center in Fig.|SJ At the edge of the masked region, 
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Fig. 6. — Reconstruction of lensing deflection angle. Left panel 
shows lensed and reconstructed unlensed CMB map; the axes are 
in arcseconds. The diagonal line shows the slice that was used to 
reconstruct the radial deflection angle shown in the right panel. 
The solid curve shows the reconstruction in the case where the 
kSZ can be perfectly removed, while the dashed line shows the 
deflection estimate when no kSZ has been removed. 



the deflection goes to zero by construction; it may be 
useful to employ an iterative approach to more accurately 
reconstruct the outer regions. The inner few arcminutcs 
are accurately reconstructed, providing information on 
the cluster mass distribution. 

A cluster will produce a second microwave background 
distortion with a near-blackbody spectrum in addition to 
lensing , namely the kinematic Suny aev-Zeldovich (kSZ) 
effect l|Sunvaev fc Zel'dovichl Il972j) . This distortion 
arises from photons scattered by the cluster gas moving 
with some bulk velocity relative to the microwave back- 
ground rest frame. If clusters were balls of gas with no 
significant internal motions, the result would be either a 
positive or negative temperature distortion proportional 
to the gas density profile, for clusters with radial ve- 
locities toward or away from the observer. This is eas- 
ily distinguishable from the dipole-like pattern induced 
by lensing. Furthermore, this spatial pattern would be 
matched by the thermal SZ signature, and could be pro- 
jected out. However, actual clusters possess significant 
internal coherent bulk velocities, particularly for highcr- 
rcdshift clusters undergoing frequent mergers, and these 
internal v elocities resu l t in m ore complex kSZ distortions 
(see, e.g.. iNaeai et alJ l|2003j) L Occasionally, the result- 
ing kSZ maps can even mimic the dipole-like pattern seen 
in lensing, although alignment with the temperature gra- 
dient would be purely coincidental. 

The degree to which the kSZ effect will degrade mass 
estimates will vary between individual clusters. Gener- 
ally, the region of significant kSZ distortion is more spa- 
tially concentrated than the lensing signal, and will have 
some correlation with the thermal SZ distortion since 
both scale with the cluster's total optical depth. Figure 
El also shows recovered deflection angles when the kSZ 
effect is included. Clearly, the lensing reconstruction is 
seriously degraded in the central cluster region where 
the kSZ signal is comparable to the lensing signal, but 
for regions further from the center, the kSZ signal does 
not alter the reconstruction. In practice, the tempera- 
ture fluctuations in an annulus from 1' to 3' away from 
the cluster will provide most of the information about 
lensing deflections. 

We hav e not included residual noise from large scale 
structure l)Vishniadll987j) . which should be at the level 
of roughly 1 /iK at the resolution of these maps. While 
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this noise makes maps less visually appealing at the low- 
est contour levels, it does not quantitatively affect the 
ability to reconstruct deflection angles. Assuming typi- 
cal gradients of 12 /zK per arcminute, noise at the level 
of 1 fiK corresponds to fluctuations at the level of 5" in 
the deflection angle. 

All of the above analysis has been for idealized maps, 
which here is equivalent to an angular resolution of 10" 
and noise levels below 1 /jK per resolution element. 
The next generation of microwave background maps at 
small scales (e.g ., Atacama Cosmology Telescope; ACT 
iKosowskvl lj2003ft or South Pole Telescope; SPT) will at- 
tain resolutions of roughly 1' and noise levels of a few /xK 
per resolution element. The ability to detect actual lens- 
ing deflections at these resolutions and sensitivities will 
be lower than presented here; quantitative analysis of 
how much the signal is degraded is ongoing. E ventually, 
the A tacama Large Millimeter Array, ALMA 3 ijWoottenl 
2002), will have both the sensitivity and resolution to 
make the kinds of maps displayed in this paper, provided 
systematic errors can be controlled well enough to allow 
reliable imaging of more than ten square arcminutes (re- 
quiring a large mosaic). 

4.3. Mass From Deflection 

In an ideal case, lensing of the microwave background 
determines the component of the deflection angle in the 
direction of the local temperature gradient. In principle, 
this is not sufficient information to determine the mass 
distribution. The deflection a(x) can be expressed as 
the gradient of a lensing potential, a — Vip; thus the 
deflection field is curl-free, V x a = 0. Thus if we in- 
fer a x from a map, in principle this determines da y /dx 
as well (though of course no practical algorithm should 
differentiate measured data!). However, we cannot ob- 
tain ct v itself by integrating with respect to because a 
different integration constant obtains for each value of y. 
So we cannot recover the complete deflection angle a nor 
the convergence (mass density) 2k — V • a, but rather 
one component of the deflection and the derivative of its 
orthogonal component. 

In practice, a number of other reasonable assumptions 
or constraints can be imposed. Assuming the deflection 
to be zero beyond some given radius will break the de- 
generacy while biasing the cluster mass estimate. The as- 
sumption of a spherically symmetric lens will also break 
the degeneracy; more generally, assuming almost any 
parametric form for the lens will largely lift the degener- 
acy. 

4.4. Errors in Cluster Catalogs 

So far we have discussed galaxy clusters individually. 
Future high-resolution microwave maps will yield sam- 
ples of hundreds to thousands of galaxy clusters selected 
by their thermal SZ spectral distortions. The lensing 
signals in these samples can be used to probe the cluster 
mass distribution in a statistical sense. In this case, it 
is important to distinguish statistical errors, which will 
average out in a large cluster sample, from systematic 
ones which will not. 

The major sources of error affecting the lensing ob- 
servations and analyses described here are residual un- 

3 www.alma.nrao.edu 



certainties from projecting out the thermal SZ spectral 
distortions and other non-blackbody foreground compo- 
nents like dust or point sources; separation of the kSZ 
signal from the lensing plus background CMB signals; er- 
rors in reconstructing the deflection map from the lensed 
CMB map; and other neglected sources of signal like the 
moving cluster effect. In order for an error to have a sys- 
tematic impact on lensing mass estimates, it must cor- 
relate with the cluster lensing signal, which is aligned 
with the local temperature gradient. This means that 
no foreground emission is likely to contribute any sys- 
tematic error to the lensing signal. The separation of 
the thermal SZ signal is also unlikely to give any sys- 
tematic error, because of its roughly symmetric shape. 
The moving cluster effect is both negligibly small for any 
single cluster and also uncorrelated with the background 
temperature gradient so it will not give any systematic 
effect over a significant survey area (i.e., a few square de- 
grees). The kinematic SZ effect can coincidentally mimic 
the lensing signal to some extent in the case of a clus- 
ter undergoing a merger, but the dipole-like morphology 
of the signal will only be aligned with the background 
temperature gradient by coincidence in a small fraction 
of these merging clusters; also any kinematic SZ dipolc- 
shaped signature will be over smaller angular scales than 
the lensing signal and will mostly be eliminated by mask- 
ing out the central portion of the cluster. 

The only clear potential systematic error arises from 
the technique used to separate the lensing signal and 
the background fluctuations. In our sample technique 
in this paper, for example, we reconstruct the deflection 
over some finite-area region around the cluster, assuming 
zero deflection at the edges of the region. This leads to 
a systematic underestimate of the deflection in the outer 
parts of the cluster, as seen in Fig-E] The extent to which 
this systematic error biases ultimate mass estimates de- 
pends on the details of the method used to reconstruct 
the deflection. This bias depends primarily on the shape 
of the background fluctuations, and for a given recon- 
struction technique can probably be modelled with a fair 
degree of certainty. Detailed estimates of the size of this 
systematic error when constructing statistical estimates 
from cluster catalogs is beyond the scope of this paper. 

5. DISCUSSION 

Cluster mass estimates based on the gas distribu- 
tion, such as those extracted from X-ray temperature 
maps or from the Sunyaev-Zeldovich effect, have in- 
trinsic uncertainties arising from the detailed relation- 
ship between the gas distribution and the total mass 
distribution, particularly at higher redshifts when clus- 
ters undergo frequent mergers. The only direct way to 
measure cluster masses is through gravitational lensing. 
Strong and weak lensing of background galaxies by a 
cluster give measurements of the shear field induced by 
the cluster mass distribution, and elaborate techniques 
have been developed to reconstruct a n estimated mass 
distri bution from lensing obse rvations i|Kaiser fc Sauiresl 
I19931 iSauires & Kaiser! Il996|) . A nu mber of cluster 
masses have been measured this way llFisher fc Ty sor 
| 1997tlKing et al.l2002UTvson et alll998HWittman et al. 
I200l . 

However, cluster mass determination via lensing of 
background galaxies has a number of limitations. At 
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usual magnitude limits, most clusters do not exhibit 
the arcs of strongly lensed background galaxies, which 
arise due to chance alignment of the cluster and back- 
ground galaxy. Mass measurements for large sam- 
ples of galaxy clusters thus cannot rely on strong lens- 
ing as a primary tool. Weak lensing of faint back- 
ground galaxies provides a more generic signal. But 
weak lensing has a number of difficult systematic con- 
siderations: (1) Orientations of background galaxies 
possess some degree o f intrinsic alignment, which can 
bias th e lensing signal llBrown et al.ll20 02l: iCatelan et all 
20011 ICritterider, et , a.U 12(1011 ICroft fc MetzleH 120011 
Heavens et al.ll2000t iPen et all 12000ft . (2) The redshift 
distribution of background galaxies is generally not well 
determined, and photometric redshi ft distributions may 
induce significant systematic errors jMcKav et al]l2002t 
ISheldon et aDl20M ISmith et al.ll200lD . C31 Thelensing 
signal becomes weaker for higher cluster redshifts be- 
cause fewer background galaxies are behind high-redshift 
clusters. (4) At a given exposure depth, weak lensing 
measurements are limited by shot noise due to a finite 
number of background galaxies. 

Lensing of the microwave background ameliorates 
these limitations. The photons all originate from a sin- 
gle, precisely-determined redshift. and this surface of 
last scattering is well behind any cluster, simplifying the 
redshift dependence of the lensing signal Furthermore, 
given sufficient resolution, the lensed temperature pat- 
tern can be probed with arbitrary accuracy, circumvent- 
ing the shot noise limits. An additional potential advan- 
tage, pointed out by Seljak and Zaldarriaga (2000), is 
that deflection angle drops off more slowly with distance 
from the cluster than does shear. However, as radius in- 
creases, the comparatively larger lensing signal is coun- 
teracted by the larger noise (arising from the primary 
CMB anisotropies) from uncertain knowledge of the un- 
lensed background; in practice this advantage over weak 
lensing is uncertain. 

The microwave background lensing signal will have its 
own systematic limitations. Foremost among these are 
the contribution of the kSZ signal and our limited knowl- 
edge about the underlying primordial temperature fluc- 
tuation pattern. As illustrated here, the impact of each 
of these varies from cluster to cluster. It is safe to say 
that some non-negligible fraction of clusters will be good 
candidates for lensing mass determinations, due to their 
location in front of a favorable area of temperature fluctu- 
ations and their lack of significant internal bulk velocities 
producing a complicated kSZ contribution. Also, which 
clusters will provide the most reliable lensing deflection 
determination will be evident from the maps themselves. 

Finally, any lensing mass estimate does not probe only 
the cluster mass, but rather all of the mass projected 
along the line of sight to the cluster. Generally, this to- 
tal projected mass is dominated by the galaxy cluster 
mass, but other masses will lead to a small systematic 
bias towards overestimating cluster masses. This can be 
addressed by comparing observations to full cosmological 
simu lations instead of to a simple cluster mass distribu- 
tion l|Chen et a l. 2003). Note that weak and strong lens- 
ing of background galaxies also suffer from this difficulty, 
although the larger redshift of the microwave background 
makes the effect somewhat more important. The linear 
contribution due to large scale structure has been crudely 



included by our use of the lensed CMB power spectrum, 
but this does not adequately include the effects of non- 
Gaussianity and non-linearity. 

Of course, the usual slew of difficulties in making ac- 
curate, high-resolution and high-sensitivity microwave 
background maps will also need to be overcome, includ- 
ing foreground emission, point sources, and the cluster's 
comparatively large thermal Sunyaev-Zeldovich distor- 
tion. Confusion-limited infrared point sources, in partic- 
ular, will be a serious problem at arcminute reso lutions 
(|BlairJll998t iBorvs et all 120031: iKnox et all 12003]) . The 
cluster thermal SZ effect can be minimized by observing 
near the frequency null around 220 GHz. Thermal SZ 
and foreground separation at the /iK level will require 
multi- frequency observations, although both will be un- 
corrected with the specific lensing temperature distor- 
tions around clusters. Accurate measurement of lensing 
signals with maximum amplitudes of a few /iK in indi- 
vidual clusters will likely require sub-arcminute resolu- 
tion observations over a range of frequencies by an in- 
strument like ALMA, and may be marginally detectable 
with ACT or SPT. 

Lensing of the CMB polarization fluctuations should 
be less susceptible to many sources of confusion than 
lensing of the temperature fluctuations considered in 
this paper: the kSZ and thermal SZ polarization fluc- 
tuations are significantly smaller than the correspond- 
ing temperature fluctuations, and infrared point sources, 
which are actually the total dust emission from high- 
redshift galaxies, are likely to show little polarized emis- 
sion. Furthermore, the primary CMB polarization pat- 
tern on small scales will have a particular spatial prop- 
erty (curl-free) that is violated in the presence of lensin g 
l|Kamionkowski et al.lll997t iZaldarriaga fc Seliaklll997[) . 
providing a "smoking gun" of the presence of lensing dis- 
tortions. If the systematic limitations to measuring lens- 
ing deflections from temperature maps prove difficult to 
overcome, polarization maps provide a viable alternative 
(although they require still greater experimental sensi- 
tivity due to the smaller polarization amplitude in the 
primary microwave background fluctuations). 

The importance of constraining galaxy cluster masses 
hardly needs to be emphasized. The number density of 
clusters of a given mass as a function of redshift is a 
highly sensitive probe of the growth rate of structure, 
and in principle could strongly constrain t he scale factor 
evolu tion at recent cosmological epochs ()Haiman et alJ 
|2001[) . The problem is that traditional cluster observa- 
tions in optical and X-ray bands, or measurements of 
the Sunyaev-Zeldovich distortion, measure the distribu- 
tion of baryons in the cluster, not its total mass. Con- 
necting the distribution of baryons and dark matter in 
clusters is a difficult and complex problem, providing 
many possibilities for systematic errors in cluster mass 
estimates, particularly at larger redshifts where clusters 
are further from dynamical equilibrium due to frequent 
merging. Lensing measurements appear to be the only 
route to a direct, reliable cluster mass estimate. A clear 
understanding and characterization of cluster lensing is 
also essential for probing the clu ster kinematic SZ signal 
l|Holderll2003l iNagai et an i2f)031. where the lensing acts 
as an unwanted complication. We hope that the images 
and arguments in this paper convincingly demonstrate 
that microwave background lensing is a viable alterna- 



tive to weak lensing of background galaxies for learning 
about the distribution of mass in galaxy clusters. 

Note: since submission of this work, two notewor- 
thy papers of re l evance have app eared in the literature 
l|Dodelsonll200llVale et alJl2004l) . These works confirm 
the results presented here and present interesting discus- 
sions of other important issues and directions for lensing 
of the CMB by galaxy clusters. 
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